The TVirbulent Heating Rate in Strong MHD Turbulence with Nonzero Cross 

Helicity 

Benjamin D. G. Chandran 

benjamin . chandran@unh . edu 
Space Science Center and Department of Physics, University of New Hampshire 

Eliot Quataert 
eliot@astro .berkeley . edu 
Department of Astronomy, University of California, Berkeley 
Gregory G. Howes 
gregory-howes@uiowa . edu 
Department of Physics & Astronomy, University of Iowa 
Joseph V. Hollweg 
joe . hollweg@unh . edu 
Space Science Center and Department of Physics, University of New Hampshire 

William Borland 
bdorland@umd . edu 
Department of Physics, University of Maryland 

ABSTRACT 

Different results for the cascade power £ in strong, incompressible MHD turbu- 
lence with nonzero cross helicity appear in the literature. In this paper, we discuss 
the conditions under which these different results are valid. Our conclusions can be 
expressed in terms of the density p, the rms amplitudes and z~ of Alfvenic fluc- 
tuations propagating parallel and anti-parallel to the background magnetic field Bq, 
and the correlation length (outer scale) measured perpendicular to Bq, denoted L^. 
We argue that if z+ ^ and if the z^ fluctuations are sustained by the reflection of 
z+ fluctuations in a strong background magnetic field, then £ ~ p(z+)^z^/L^ as in 
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previous studies by Hossain, Matthaeus, Dmitruk, Lithwick, Goldreich, Sridhar, and 
others. On the other hand, if the minority wave type (z ) is sustained by some form 
of forcing that is uncorrelated with or only weakly correlated with the z+ fluctuations, 
then £ can be much less than p(z+)^z^/L^, as in previous studies by Dobrowolny, 
Lazarian, Chandran and others. The mechanism for generating the minority wave type 
strongly affects the cascade power because it controls the coherence time for interac- 
tions between oppositely directed wave packets at the outer scale. 

Subject headings: turbulence — magnetic fields — magnetohydrodynamics — solar 
wind — solar corona — solar flares 



Turbulence plays an important role in numerous astrophysical plasmas, from the solar corona 
and solar wind (Coleman 1968; Belcher & Davis 1971; Tu & Marsch 1995; Goldstein, Matthaeus, 
& Roberts 1995) to star-forming molecular clouds and clusters of galaxies (Boldyrev, Nordlund, 
& Padoan 2002; Narayan & Medvedev 2001; Cho & Lazarian 2003; Elmegreen & Scalo 2004). 
In general, astrophysical turbulence is compressible, and the compressive component of the turbu- 
lence may play an important role in processes such as ion heating and the scattering of energetic 
particles (Li & Habbal 2001; Yan & Lazarian 2004; Chandran 2005, 2008a). However, in this pa- 
per we focus on the non-compressive component of the turbulence, which generally dominates the 
energy at scales smaller than the coUisional mean free path, since (non-compressive) Alfven waves 
experience negligible linear coUisionless damping at wavelengths X much larger than the proton 
gyroradius and frequencies (O much smaller than the proton cyclotron frequency Q, (Barnes 
1966). 

For X ^ n and co <^ Q.u the non-compressive component of the turbulence can be modeled 
using incompressible magnetohydrodynamics (MHD), even in a low-coUisionality plasma. [See 
Howes et al (2008) and Schekochihin et al (2009) for recent discussions of this point.] In incom- 
pressible MHD, the fluctuating velocity and magnetic fields can be decomposed into fluctuations 
propagating parallel and anti-parallel to the background magnetic field Bq, denoted z+ and re- 
spectively. When the energies in z+ and fluctuations differ, the turbulence has nonzero cross 
helicit\n and is referred to as "imbalanced." 



Turbulence in the solar corona and inner solar wind is imbalanced, with an excess of Alfvenic 
fluctuations propagating away from the Sun (Roberts et al 1987; Grappin et al 1990; Bavassano 



The cross helicity is defined as / d\ v • B, where v is the velocity and B is the magnetic field. 



1. Introduction 
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et al 2000a, 2000b; Bruno & Carbone 2005). In order to model the heating of the solar corona and 
solar wind, it is of interest to determine the turbulent heating rate or cascade power in imbalanced 
MHD turbulence. Different expressions for this cascade power have appeared in the literature (see, 
e.g., Dobrowolny, Mangeney, & Veltri 1980; Hossain et al 1995; Lithwick, Goldreich, & Sridhar 
2007; Beresnyak & Lazarian 2008; Chandran 2008b). Our goal in this paper is to clarify the 
conditions under which these different expressions are valid. We focus on the difference between 
situations in which the minority Alfven wave type (which we take to be z^) is sustained by the 
reflection of z+ fluctuations and situations in which the z fluctuations are sustained by some form 
of forcing. Although a number of studies have addressed the effects of cross helicity on the inertial- 
range power spectra of z+ and fluctuations (e.g., Grappin, Pouquet, & Leorat 1983; Lithwick 
et al 2007; Beresnyak & Lazarian 2008, 2009; Chandran 2008b; Perez & Boldyrev 2009; Podesta 
& Bhattacharjee 2009), we do not address the inertial-range power spectra in this paper. 



2. Theoretical Framework 

We consider Alfvenic (i.e., transverse and non-compressive) fluctuations in a stationary, inho- 
mogeneous medium, with background density p and background magnetic and velocity fields Bq 
and u = uBq/Bq. The velocity and magnetic-field fluctuations 6v and 5B can be written in terms 
of the Elsasser variables 

z± = 5vT^. (1) 

The equations for z± are (Velli 1993, Verdini & Velli 2007) 

^ + (u±va)-Vz±+zT-V(utva)t^(z±-zT) (^V-vaT^V-u^ =-zT-Vz±-^, 

(2) 

where va = Bo/y/Anp and jCtot = p+B^ /Sn is the sum of the thermal and magnetic pressure. The 
second term on the left-hand side of equation (O represents wave propagation at velocity ±va 
in the frame of the plasma. The third and fourth terms on the left-hand side give rise to wave 
reflection. The right-hand side of equation ^ contains the nonlinear terms. Nonlinear interactions 
arise only when oppositely propagating wave packets interact (Iroshnikov 1963; Kraichnan 1965). 
In this paper, we focus on turbulent fluctuations at the "outer scale" or energy-injection scale. We 
view these outer- scale fluctuations as a collection of z+ and z^ wave packets with correlations 
lengths and Ljj^ perpendicular and parallel to Bq, where for simplicity we have taken to be 
the same for both the z+ and z fluctuations. 

The shearing of a z^ wave packet at the outer scale is dominated by interactions with z^ 
wave packets at the outer scale. Since we take the fluctuations to be Alfvenic in nature, z^ is 
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perpendicular to Bq, and the magnitude of the ■ Vz^ temi in equation ([21) is ~ ZqZq /L^, where 
Zq is the rms amplitude of the fluctuations at the outer scale. The rate of shearing of z"^ wave 
packets by z^ wave packets at the outer scale is thus 



■'shear ~ r " 



Nonlinear interactions between z+ and z^ fluctuations can be thought of as "collisions" between 
oppositely directed wave packets. If we consider some point P moving with a wave packet at 
velocity ±vaz, the time required for this point to pass through an outer-scale z^ wave packet is 
approximately 

a 

^coll ~ 



^.^oi, ~ — (4) 



VA 

(where we have ignored the factor of 1/2 that results from the fact that the relative speed of point P 
and the moving z^ wave packet is 2va). We define the quantity 

ZqL^ 

L^Va 

If %^ < 1, then %^ is approximately the fractional change in an outer-scale wave packet caused 
by shearing by outer-scale fluctuations during a single "collision" of duration If %^ > 1, 
then an outer- scale z^ wave packet undergoes a fractional change of order unity in a time ~ L^/z^, 
which is < t^^^^. 

We take z^ to be the dominant wave type, and assume that 

Zq<^z^<va. (6) 

which implies that 55 < Bq. Throughout this paper, we focus on the strong-turbulence limit, in 
which 

X+>1- (7) 
We note that the conditions %^ > 1 and 55 -C Bq are compatible if L^/Lt is sufficiently small. 



3. Forced Turbulence 

If a magnetized plasma is physically stirred by a set of objects distributed throughout its vol- 
ume, and if the motions of the different stirring objects are uncorrected, then this stirring will 
excite both z^ and z^ waves, with comparable power going into each wave type. If this stirring 
is the primary mechanism for exciting Alfvenic fluctuations, then the resulting turbulence will be 
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approximately balanced, with little cross helicity. However, volumetric forcing can play an impor- 
tant role in imbalanced turbulence if, in addition to the forcing, a large flux of the dominant z+ 
fluctuations is launched into the plasma from the plasma boundary. If the forcing is weak, then Zq 
will be small, the cascade time of the z+ fluctuations will be large, and the z+ fluctuations injected 
at the plasma boundary will propagate a large distance into the plasma before their energy cas- 
cades. The role of forcing is then different for z+ and ■ Within a certain distance of the boundary 
(see section 13.11) . forcing makes only a small contribution to the energy of the z+ fluctuations, 
but is the mechanism that sustains the z fluctuations. We note that if z+ fluctuations satisfying 
equation ^ are launched into the plasma from the plasma boundary, and if z energy is not con- 
tinually injected into the system throughout its volume, then z^ energy will decay on the short 
time scale L^/zq, leaving the system in a maximally aligned state|^ with vanishing cascade power 
(Dobrowolny et al 1980; Matthaeus & Montgomery 1980; Pouquet, Sulem, & Meneguzzi 1988). 

We model the forcing of z^ fluctuations as a stochastic source term F(x, t) added to right-hand 
side of equation For simplicity, we take the spatial variation of F in directions perpendicular 
to Bq to be characterized by a single scale length (correlation length), which we set equal to the 
perpendicular correlation length of the outer-scale z+ fluctuations, L±. We also take the spatial 
variation of F along Bq to be characterized by a single scale length, denoted L|| f, which can differ 
from L|j". We assume that the (stochastic) time variation of F(x,?) is characterized by the single 
time scale L|| f/vA- 

The properties of outer-scale z^ wave packets are determined by both (1) the way in which 
the z^ fluctuations are sheared by z+ fluctuations and (2) the source term. We discuss these two 
mechanisms in turn and their consequences for the correlation time of the shearing that is applied 
by z^ wave packets and experienced by a single, moving z+ wave packet. Because of equation (|7]), 
z^ fluctuations can only propagate a distance < L'^ before their energy cascades to smaller scales. 
As a result, if we consider the outer-scale z+ wave packet labeled "A" in Figure [H the outer-scale 
z^ fluctuations encountered by wave packet A are sheared primarily by wave packet A and the 
outer-scale wave packet immediately "ahead" of wave packet A, which is labeled wave packet B 
in Figure [B Since wave packet B remains in front of wave packet A as the two wave packets 
propagate along the magnetic field, the shearing experienced by the outer-scale z^ fluctuations 
that are encountered by wave packet A remains coherent throughout the cascade time of wave 
packet A. On the other hand, the properties of the source term F in the region that contributes to 
the z" fluctuations encountered by wave packet A decorrelates over a time scale ~ L|| f/vA, because 
of both the space and time dependence of F. Thus, the coherence time for the shearing experienced 



^When z = 5v + SB/y/4-np — 0, the state is called "maximally aligned" because 5v and 5B are everywhere anti- 
parallel. 
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by the outer-scale z+ wave packets 



■^cor 



VA 



(8) 



We note that if L|| f ^ L^, then because of equation (|7]) the shearing experienced by outer-scale 
fluctuations will cause the fluctuations to vary strongly over a distance along the field, 
causing to become comparable to L^. However, in this case, the shearing experienced by the z+ 

fluctuations remains coherent for a time L|| f/vA rather than the collision time scale ~ /vp^, 
for the reasons given above. 



B 



a 



-ii,f 



Fig. 1 . — A and B represent two outer-scale z+ wave packets propagating along the background magnetic field Bq. 
In this pictorial example, the parallel correlation length of the forcing term, L|| f is taken to be longer than the parallel 
correlation length of the outer-scale z+ wave packets, Li . 



We define the quantity 



We first consider the case in which 



1 



vaL± ' 



(9) 



r < 1- (10) 

When equation (flOl) is satisfied, x is approximately the fractional change in an outer-scale z+ 
wave packet during a time t^^^^. During successive time intervals of duration t^T^^^, the incremental 
changes in the z+ wave packet that are caused by the shearing of fluctuations add like a random 
walk. Thus, (%^) time intervals of duration r^^j. are required in order for the incremental changes 
in the z+ wave packet to accumulate to an order-unity fractional change. The cascade time for 
outer- scale z+ wave packets is then given by 



L2 



■^casc 



VA 



(11) 



^Equation ^ is valid unless L|| f/vA exceeds the z+ cascade time f, 
time and cascade power are given by equations ( fT4] i. ( fTSI l, and (TH 



case in which case f^oi- 



' f the cascade 
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e+ ~ • ^""^ ^""^ '"^ ■ (12) 



The cascade power in z+ fluctuations, £+ ~ Pl^o )^/^casc' then (Beresnyak & Lazarian 2008, 
Chandran 2008b) 

vaL]_ 

If we set L± = L[| f, then equation (fT2)) reduces to the resuh of Dobrowolny et al (1980) for 
isotropic turbulence. In the case that L|| f ~ L|j" and ~ 1, we can rewrite equation (fT2l) as 
(Chandran 2008b) 

e+^PiSl!. (13) 

Because of equation ([7]), the cascade time for fluctuations is simply t^^^^. ~ L^/zq, and the 
cascade power in z^ fluctuations is 

^ p4(^o)^ 



La 

comparable to the value of £+ in equation (fT3l) .R 



(14) 



If, instead of equation (ITOI) . 

r > 1 (15) 

then the shearing experienced by outer-scale z+ wave packets remains coherent throughout their 
cascade time, 

~ — (16) 

^0 

and 

8+ ~ P^^l^, (17) 
L± 

while £ is still given by equation (fT4l) . Turbulent heating rates of this form were previously ob- 
tained by Hossain et al (1995) on the basis of numerical simulations of decaying MHD turbulence, 
and by Lithwick et al (2007) on the basis of theoretical arguments. 



^In the weak-turbulence limit X ^ 1' both e+ and scale in the same way with Zq, Lj_, and L|| [in 

particular, 8* °^ p{zq )^{zq )^L\\/{vaL]_) (Lithwick & Goldreich 2003)]. However, there is an additional dependence 
of on the inertial -range spectral indices, so that the fluctuation type (z+ or z^) with the steeper inertial -range 
power spectrum has the larger cascade power. It was shown analytically by Galtier et al (2000) and Lithwick & 
Goldreich (2003) that this dependence becomes very strong as the spectral indices approach the limiting values at 
which the assumption of local interactions breaks down, so that the ratio e+/e^ can take on any value. In the strong 
turbulence limit X~ ^ ~ 1' ^"d also scale in the same way with Zq and Lj_ when L|| f ^ and x < 
1, as shown in equations ( fT3] l and ( fT4] l. However, if e+ and depend on the inertial-range spectral indices in 
the strong turbulence case as they do in the weak-turbulence limit, then additional coefficients (not necessarily of 
order unity) need to be included in equations (T3[ and (fT4l i. The presence of such coefficients would allow for the 
possibility of forced, steady-state, strongly turbulent flows in which e+/e^ differs significantly from unity, even when 
the assumptions leading to equations (13[ and (fT4l i are valid (i.e., x^ — 1, X < 1' ™d Lii f ^ Lf). 
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3.1. An Example of Forced, Imbalanced MHD Turbulence 

In this section, we consider a finite plasma that is stirred by a set of objects located throughout 
its volume. We take the dimension and spacing of these objects perpendicular to the magnetic field 
to be L^, and we take the objects to be spaced by a distance along the magnetic field. We take 
the correlation time of the forcing to be L|j"/vA. In the absence of wave launching at the boundary, 
these stirring objects set up balanced MHD turbulence, in which the rms amplitudes of the z+ and 
z fluctuations are both equal to a quantity that we denote 5z. We define 

10 = (18) 

L^VA 

and assume that 

Xo«l- (19) 
Thus, in the absence of wave launching, the forcing would set up weak turbulence at the outer 
scale, with an outer-scale cascade time of (%o)^^^|[^/va and a cascade power of 

p(5z)'U 

(20 ~ -72 -■ (20) 

L^va 

We now assume that z+ waves are launched into the plasma from one its boundaries (the 
"launching boundary") with correlation lengths and perpendicular and parallel to the mag- 
netic field. We consider the region that is sufficiently close to the launching boundary that the 
amplitudes of the z+ fluctuations have not decayed very much from their values right next to the 
boundary. (The size of this region is estimated below.) We assume that within this region, 

X+-1, (21) 

implying that Zq ^ 5z. We also assume that the presence of the additional z+ waves does not sub- 
stantially alter the rate at which energy is injected into either the z+ of fluctuations!^ Since the 
cascade time for the z~ fluctuations is short, the rate at which energy is injected into z fluctuations 
will be in balance with the rate at which z^ energy cascades to small scales. Thus, 

Qq ~ — • (^^) 



^This assumption is plausible for the following reason. The net power per unit volume going into the fluctu- 
ations, denoted Q^, is approximately (F • z*), where (. . .) represents an ensemble average. The value of can be 
thought of as consisting of two components, z* = z]*^ + Z2 , where is correlated with F, and Zj is not. A rough 
estimate of is then frmszfrms' where Frms and zf^n^^ are the rms amplitudes of F and . Because the correlation 
time of F in our example is LJ^/va, which is comparable to or shorter than the cascade time of each fluctuation type, 
rms ■fimsi'^ / for both wave types, so that ^ Q . 
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Equations (|20l) through (|22l) then lead to the relation 



+ , (23) 



which illustrates how the enhanced level of z+ fluctuations reduces Zq relative to its value (6z) in 
the absence of wave launching from the boundary. 



Determining the rate at which z energy decays in this example is subtle. Equations (|T2l) and 
(1231) imply that 

£+ ~ —2 (24) 

L^VA 

which is comparable to the rate at which energy is injected by the forcing term. The question of 
whether the z+ energy increases or decreases with increasing distance from the launching boundary 
thus depends upon the undetermined factors of order unity in equations (|20l) and ((24l) . To proceed, 
we conjecture that the z+ and power spectra are "pinned" (i.e., equal) at the dissipation scale 
(Grappin et al 1983; Lithwick & Goldreich 2003 — but see Perez & Boldyrev 2009 and Beresnyak 
& Lazarian 2009), and that there is a dimensionless prefactor in the expression for £+ (or £^) 
that increases as the z+ (or z^) inertial-range power spectrum becomes steeper, as in the case of 
imbalanced weak MHD turbulence (Galtier et al 2000; Lithwick & Goldreich 2003). This implies 
that in the imbalanced case in which Zq ^ Zq, the value of £+ will be somewhat larger than both £^ 
and Qq, so that the z+ energy will decay in time, when viewed from a frame of reference moving 
with the z+ fluctuations at velocity vaz. The time scale r+j,^ for this decay is given approximately 
by equation (fTTl) . with L|| f set equal to LjJ^. In the plasma frame, the z+ energy then decays over a 
distance of L+^^y ~ va^+sc or (since %+ ~ 1) 



^d+ecay-^Xo'- (25) 



At distances < L^^^^y from the launching boundary, the turbulence would then be highly imbal- 
anced, with forcing sustaining the z fluctuations but making only a small contribution to the 
energy in z+ fluctuations. On the other hand, at distances ^ ^jecay from the launching boundary 
the turbulence would become approximately balanced. If our conjecture above about the relation 
between the cascade power and inertial-range spectral index in strong imbalanced turbulence is 
incorrect, then other phenomenologies are possible. 



4. Reflection-Driven Turbulence in a Strong Background Magnetic Field 



In this section, we assume that the z wave packets are generated by reflection of the dominant 
z+ wave packets. The key difference from the last section is that the mechanism generating the 
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z fluctuations is now strongly correlated with, and in fact determined by, the properties of the z+ 
fluctuations. We begin by considering outer-scale z+ fluctuations that satisfy the "critical-balance 
condition" (Goldreich & Sridhar 1995; Lithwick et al 2007), 

X+~l. (26) 

We consider the case 1 at the end of this section. Part of our motivation for considering 

outer scales satisfying equation (|26|) comes from the work of Cranmer & van Ballegooijen (2005), 
who modeled the launching of Alfven waves into the corona and solar wind from the Sun. They 
based their analysis on observations of photospheric motions and a model of wave propagation in 
the chromosphere, corona, and solar wind. They found that for heliocentric distances between 2 
and 10 solar radii, the bulk of the wave power is at periods of 1 — 5 minutes, with wave amplitudes 
such that equation (l26l) is approximately satisfied. 

Equation (|26l ) implies that the outer-scale wave packets only travel a distance ~ L|j^ before 
their energy cascades to smaller scales. Thus, the outer-scale fluctuations present at an arbitrary 
point P originated within a distance ~ L|j^ from P. If we consider two outer-scale z~^ wave packets 
A and B of length L|j^ and width L±, with A right behind B as in Figure [21 then the outer-scale 
z^ wave packets encountered by wave-packet A (denoted B' in Figure |2l) will primarily be the z^ 
field arising from the reflection of wave packet B. If we assume that the gradient scale lengths of 
the background medium are at least as long as the "damping length" VA^casc of the z+ wave packet, 
then the z^ wave packets encountered by wave packet A will continue to have the same form and 
structure until wave packet B changes. The shearing experienced by wave packet A is thus coherent 
throughout its cascade time, rather than a random-walk-like accumulation of uncorrelated, small- 
amplitude distortions. The cascade time of wave packet A is therefore simply ~ ((Ogjjg^.)^^ or 

~ — (27) 

'case — ' \^ ' J 

^0 

which implies that 

e ~ — , (28) 

as in the studies of Hossain et al (1995) and Lithwick et al (2007). 

Lithwick et al (2007) argued that in strong, imbalanced, MHD turbulence, the shearing felt by 
the z+ wave packets throughout the inertial range is coherent over a time that significantly exceeds 
the wave-packet "collision time" or crossing time, independent of the mechanism that generates 
the outer-scale fluctuations. Our arguments differ from theirs in two ways. First, our discussion 
is restricted to the outer-scale fluctuations. Second, we argue that equation (l27l) is valid when the 
minority z^ wave type is generated by wave reflection, but not when z^ is generated by forcing 
that is only weakly correlated with the z+ fluctuations (unless > 1, as described in section[3]). 
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Fig. 2. — In strong, reflection-driven turbulence, the properties of the wave packet (B') encountered by the z+ 
wave packet on the left (A) are determined by the wave packet on the right (B). The shearing of wave packet A by 
wave packet B' thus remains coherent over a time comparable to the cascade time of wave packet B. 

For the waves launched into the solar corona by photospheric motions, it is typically argued 
that ~ 10^ cm at the coronal base, and that the Alfven waves at this perpendicular scale have 
a broad frequency spectrum, extending up to wave periods of many hours (Dmitruk & Matthaeus 
2003, Cranmer & van Ballegooijen 2005, Verdini & Velli 2007). Although Cranmer & van Balle- 
gooijen (2005) find that the energy in outward Alfven waves (z+ in our discussion) in the corona 
is dominated by waves with periods of ~ 1 — 5 minutes, wave reflection is more efficient for low- 
frequency waves. Thus, if nonlinear wave- wave interactions were neglected, then the inward wave 
packets would have a longer parallel correlation length than the outward waves. 

We now explore the effects of frequency-dependent reflection, including the effects of nonlin- 
ear wave-wave interactions, by considering the limiting case in which the source term for inward 
{z~) waves is invariant along the magnetic field. The parallel correlation length of the fluctu- 
ations is given by the arguments of Lithwick et al (2007), as follows. We continue to assume the 
validity of equation (|26l) . which is motivated for the coronal case by the results of Cranmer & van 
Ballegooijen's (2005) model of wave launching and propagation. Equation (|26|) implies that a 
wave packet can only travel a distance of order before the wave packet is strongly distorted. 
Since the z+ fluctuations decorrelate over a distance along the magnetic field, the inward waves 
acquire a parallel correlation length of by interacting with the outward waves, even though the 
parallel correlation length of the source term for the inward waves is ^ L|j^. 

Nevertheless, the coherence time for the shearing experienced by an outer-scale z+ wave 
packet remains much longer than This can again be seen with the aid of Figured Now, the 
source term for the inward waves is invariant along the magnetic field, but the properties of the 
wave packet (B') encountered by the z"*" wave packet A are determined by the way in which 
the z+ wave packet B shears the inward waves generated by the uniform source term. Since wave 
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packet A and B remain adjacent as they propagate outward along the magnetic field, the z~ wave 
field encountered by wave packet A will remain approximately the same until wave packet B 
changes, i.e., for a time equal to the cascade time of the outward waves. Since the shearing of the 
z+ wave packets is coherent throughout their cascade time, their cascade time is again given by 
equation (ITTI) . and the cascade power for z+ fluctuations is again given by equation (l28l) . 

We now turn briefly to the case in which ^ 1 , with the fluctuations generated by wave 
reflection. As in the case ~ 1, the shearing experienced by outer-scale z+ wave packets remains 
coherent throughout their cascade time, and the cascade power in z+ fluctuations is again given by 
equation dlH]), as in the studies by Matthaeus et al (1999) and Dmitruk & Matthaeus (2003). 

Although we have focused on wave reflection in this section, similar arguments would apply 
to a compressible plasma in which the minority waves (z~) are generated by the parametric decay 
of z+ waves (Galeev & Oraevskii 1963; Sagdeev & Galeev 1969; Vinas & Goldstein 1991; Malara 
et al 2000, 2001; Del Zanna et al 2001). In this case, the waves encountered by the z+ wave 
packet A in Figure [2] will again be determined primarily by the properties of the z+ wave packet B, 
which is immediately ahead of wave packet A, and the shearing experienced by wave packet A 
will remain coherent throughout the cascade time of wave packet A. 



5. Application to the Solar Corona and Solar Wind 

In the solar corona, the dominant source of low-frequency. Sunward-propagating, Alfven 
waves is likely wave reflection (Matthaeus et al 1999; Dmitruk et al 2001, 2002). The anti-Sunward 
waves launched by photospheric motions have sufficiently large amplitudes and low frequencies 
that X^ ^ 1 within the corona (Cranmer & van Ballegooijen 2005; Cranmer, van Ballegooijen, & 
Edgar 2007; Verdini & Velli 2007). If these anti-Sunward waves have enough time to develop a 
turbulent power spectrum extending to the (perpendicular) dissipation scale before they propagate 
out of the corona to larger heliocentric distances, then the cascade power and turbulent heating rate 
in the corona are given by equation (|28] ). as in the studies of Matthaeus et al (1999), Dmitruk & 
Matthaeus (2003), Cranmer & van Ballegooijen (2005), and Verdini & VeUi (2007). The appUca- 
tion of equations (fT2l) and (fT3] ) to the corona by Chandran (2008b) was incorrect. The use of the 
correct equation in models of the coronal is critical, because the energy of anti-Sunward waves dra- 
matically exceeds the energy of Sunward waves in the corona (Cranmer & van Ballegooijen 2005; 
Verdini & Velli 2007). When the correct expressions for e+ and are used in models of the 
turbulent heating of open-magnetic-field-line regions in the solar corona, the heating rate is ener- 
getically sufficient to explain the acceleration of the fast solar wind for plausible assumptions about 
the amplitude and frequency spectrum of the waves launched from the Sun (Matthaeus et al 1999; 
Dmitruk et al 2002; Cranmer & van Ballegooijen 2005; Verdini & Velli 2007). On the other hand. 
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if equations (fT2l) and/or (fT3]) were used, then the heating rate in these models would be much too 
small to explain the acceleration of the fast wind, because Zq /zq is so small in the corona. [See 
Cranmer & van Ballegooijen (2005) and Verdini & Velli (2007).] 

In the solar wind at low heliographic latitudes, the Sunward-propagating Alfven waves may be 
driven primarily by instabilities resulting from shear in the solar-wind velocity (Roberts et al 1992; 
Breech et al 2008). It is not entirely clear how to model this process, but a reasonable first ap- 
proximation is to treat the generation of Sunward waves by velocity-shear instabilities as a form 
of mechanical forcing that is only weakly correlated with the anti-Sunward waves. If we denote 
the parallel correlation length of this forcing as L|| f, and if %~ < 1, then the cascade power and 
turbulent heating rate are given by equation (fT2l) . which is a factor of % smaller than the value 
in equation ([28]) (Beresnyak & Lazarian 2008; Chandran 2008b). If L|| f ~ L|[" and ~ 1, then 
the cascade power and turbulent heating rate are given by equation (fT3l) (Chandran 2008b). The 
difference between the different expressions for the turbulent heating rate in reflection-driven tur- 
bulence and forced turbulence makes less of a difference in models of the solar wind at heliocen- 
tric distances of 1 AU and beyond than in models of the corona, because Zq /zq is not so small. 
Nevertheless, this difference may be significant when comparing model temperature profiles with 
observations, as shown in the recent study by Ng et al (2009). 



6. Conclusion 

In this paper, we have argued that the cascade power in strong, imbalanced, incompressible 
MHD turbulence depends upon the mechanism that sustains the minority Alfven wave type (z^) 
at the outer scale. If the outer-scale z^ wave packets are sustained by the reflection of z+ wave 
packets, then the shearing of the outer-scale z+ wave packets remains coherent throughout their 
cascade time, a time that can greatly exceed the crossing time of two counter-propagating outer- 
scale wave packets. On the other hand, if the outer-scale wave packets are sustained by some 
form of forcing that is only weakly correlated with the z+ fluctuations, then the coherence time for 
the shearing of the outer-scale z+ wave packets is shorter than the cascade time, unless the parallel 
correlation length of the forcing function is very long (i.e., unless % > 1). The longer coherence 
time in reflection-driven turbulence enhances the cascade power in reflection-driven turbulence 
relative to forced turbulence with the same values of Zq,L^, and L^. 
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